FERMILAB-Pub-98/091-E 
Phys. Rev. D 58, 072001 (1998) 

Observation of B + -> ip(2S) K+ and B° -> ^(25) K*(892)° decays 
and measurements of B-meson branching fractions into J/ift and 

ift(2S) final states 



66 "■ 

G\ ; F. Abe, 17 H. Akimoto, 39 A. Akopian, 31 M. G. Albrow, 7 A. Amadon, 5 S. R. Amendolia, 27 

^ ; D. Amidei, 20 J. Antos, 33 S. Aota, 37 G. Apollinari, 31 T. Arisawa, 39 T. Asakawa, 37 

^ W. Ashmanskas, 18 M. Atac, 7 P. Azzi-Bacchetta, 25 N. Bacchetta, 25 S. Bagdasarov, 31 

M. W. Bailey, 22 P. de Barbara, 30 A. Barbaro-Galtieri, 18 V. E. Barnes, 29 B. A. Barnett, 15 
! M. Barone, 9 G. Bauer, 19 T. Baumann, 11 F. Bedeschi, 27 S. Behrends, 3 S. Belforte, 27 

Q\ [ G. Bellettini, 27 J. Bellinger, 40 D. Benjamin, 35 J. Bensinger, 3 A. Beretvas, 7 J. P. Berge, 7 

J. Berryhill, 5 S. Bertolucci, 9 S. Bettelli, 27 B. Bevensee, 26 A. Bhatti, 31 K. Biery, 7 
^ ■ C. Bigongiari, 27 M. Binkley, 7 D. Bisello, 25 R. E. Blair, 1 C. Blocker, 3 S. Blusk, 30 



J. Berryhill, 5 S. Bertolucci, 9 S. Bettelli, 27 B. Bevensee, 26 A. Bhatti, 31 K. Biery, 7 
C. Bigongiari, 27 M. Binkley, 7 D. Bisello, 25 R. E. Blair, 1 C. Blocker, 3 S. Blusk, 30 
. Bodek, 30 W. Bokhari, 26 G. Bolla, 29 Y. Bonushkin, 4 D. Bortoletto, 29 J. Boudreau 
Breccia, 2 C. Bromberg, 21 N. Bruner, 22 R. Brunetti, 2 E. Buckley-Geer, 7 H. S. Budd, 30 
3urkett, 20 G. Busetto, 25 A. Byon- Wagner, 7 K. L. Byrum, 1 M. Campbell, 20 A. Caner, 27 
W. Carithers, 18 D. Carlsmith, 40 J. Cassada, 30 A. Castro, 25 D. Cauz, 36 A. Cerri, 27 
. Chang, 33 P. T. Chang, 33 H. Y. Chao, 33 J. Chapman, 20 M. -T. Cheng, 33 M. Chertok, 34 
Chiarelli, 27 C. N. Chiou, 33 F. Chlebana, 7 L. Christofek, 13 M. L. Chu, 33 S. Cihangir, 7 

A. G. Clark, 10 M. Cobal, 27 E. Cocca, 27 M. Contreras, 5 J. Conway, 32 J. Cooper, 7 
ivi. Cordelli, 9 D. Costanzo, 27 C. Couyoumtzelis, 10 D. Cronin-Hennessy, 6 R. Culbertson, 5 

D. Dagenhart, 38 T. Daniels, 19 F. DeJongh, 7 S. DelPAgnello, 9 M. Dell'Orso, 27 R. Demina, 7 
L. Demortier, 31 M. Deninno, 2 P. F. Derwent, 7 T. Devlin, 32 J. R. Dittmann, 6 S. Donati, 27 
' Done, 34 T. Dorigo, 25 N. Eddy, 20 K. Einsweiler, 18 J. E. Elias, 7 R. Ely, 18 E. Engels, Jr., 28 

T. Erdmann, 7 D. Errede, 13 S. Errede, 13 Q. Fan, 30 R. G. Feild, 41 Z. Feng, 15 C. Ferretti, 27 

I. Fiori, 2 B. Flaugher, 7 G. W. Foster, 7 M. Franklin, 11 J. Freeman, 7 J. Friedman, 19 
Fukui, 17 S. Gadomski, 14 S. Galeotti, 27 M. Gallinaro, 26 O. Ganel, 35 M. Garcia-Sciveres, 18 
A. F. Garfinkel, 29 C. Gay, 41 S. Geer, 7 D. W. Gerdes, 15 P. Giannetti, 27 N. Giokaris, 31 
P. Giromini, 9 G. Giusti, 27 M. Gold, 22 A. Gordon, 11 A. T. Goshaw, 6 Y. Gotra, 25 
K. Goulianos, 31 H. Grassmann, 36 L. Groer, 32 C. Grosso-Pilcher, 5 G. Guillian, 20 
J. Guimaraes da Costa, 15 R. S. Guo, 33 C. Haber, 18 E. Hafen, 19 S. R. Hahn, 7 R. Hamilton, 11 
T. Handa, 12 R. Handler, 40 F. Happacher, 9 K. Hara, 37 A. D. Hardman, 29 R. M. Harris, 7 

F. Hartmann, 16 J. Hauser, 4 E. Hayashi, 37 J. Heinrich, 26 W. Hao, 35 B. Hinrichsen, 14 
K. D. Hoffman, 29 M. Hohlmann, 5 C. Hoick, 26 R. Hollebeek, 26 L. Holloway, 13 Z. Huang, 20 

B. T. Huffman, 28 R. Hughes, 23 J. Huston, 21 J. Huth, 11 H. Ikeda, 37 M. Incagli, 27 
J. Incandela, 7 G. Introzzi, 27 J. Iwai, 39 Y. Iwata, 12 E. James, 20 H. Jensen, 7 U. Joshi, 7 

E. Kajfasz, 25 H. Kambara, 10 T. Kamon, 34 T. Kaneko, 37 K. Karr, 38 H. Kasha, 41 Y. Kato, 24 
T. A. Keaffaber, 29 K. Kelley, 19 R. D. Kennedy, 7 R. Kephart, 7 D. Kestenbaum, 11 

D. Khazins, 6 T. Kikuchi, 37 B. J. Kim, 27 H. S. Kim, 14 S. H. Kim, 37 Y. K. Kim, 18 L. Kirsch, 3 
S. Klimenko, 8 D. Knoblauch, 16 P. Koehn, 23 A. Kongeter, 16 K. Kondo, 37 J. Konigsberg, 8 
K. Kordas, 14 A. Korytov, 8 E. Kovacs, 1 W. Kowald, 6 J. Kroll, 26 M. Kruse, 30 



W 
] 

A 



53 V 



1 



S. E. Kuhlmann, 1 E. Kuns, 32 K. Kurino, 12 T. Kuwabara, 37 A. T. Laasanen, 29 I. Nakano, 12 

S. Lami, 27 S. Lammel, 7 J. I. Lamoureux, 3 M. Lancaster, 18 M. Lanzoni, 27 G. Latino, 27 
T. LeCompte, 1 S. Leone, 27 J. D. Lewis, 7 P. Limon, 7 M. Lindgren, 4 T. M. Liss, 13 J. B. Liu, 30 
Y. C. Liu, 33 N. Lockyer, 26 O. Long, 26 C. Loomis, 32 M. Loreti, 25 D. Lucchesi, 27 P. Lukens, 7 

S. Lusin, 40 J. Lys, 18 K. Maeshima, 7 P. Maksimovic, 19 M. Mangano, 27 M. Mariotti, 25 
J. P. Marriner, 7 A. Martin, 41 J. A. J. Matthews, 22 P. Mazzanti, 2 P. Mclntyre, 34 P. Melese, 31 
M. Menguzzato, 25 A. Menzione, 27 E. Meschi, 27 S. Metzler, 26 C. Miao, 20 T. Miao, 7 
G. Michail, 11 R. Miller, 21 H. Minato, 37 S. Miscetti, 9 M. Mishina, 17 S. Miyashita, 37 
N. Moggi, 27 E. Moore, 22 Y. Morita, 17 A. Mukherjee, 7 T. Muller, 16 P. Murat, 27 S. Murgia, 21 
H. Nakada, 37 I. Nakano, 12 C. Nelson, 7 D. Neuberger, 16 C. Newman-Holmes, 7 
C.-Y. P. Ngan, 19 L. Nodulman, 1 A. Nomerotski, 8 S. H. Oh, 6 T. Ohmoto, 12 T. Ohsugi, 12 

R. Oishi, 37 M. Okabe, 37 T. Okusawa, 24 J. Olsen, 40 C. Pagliarone, 27 R. Paoletti, 27 
V. Papadimitriou, 35 S. P. Pappas, 41 N. Parashar, 27 A. Parri, 9 J. Patrick, 7 G. Pauletta, 36 
M. Paulini, 18 A. Perazzo, 27 L. Pescara, 25 M. D. Peters, 18 T. J. Phillips, 6 G. Piacentino, 27 

M. Pillai, 30 K. T. Pitts, 7 R. Plunkett, 7 A. Pompos, 29 L. Pondrom, 40 J. Proudfoot, 1 
F. Ptohos, 11 G. Punzi, 27 K. Ragan, 14 D. Reher, 18 M. Reischl, 16 A. Ribon, 25 F. Rimondi, 2 
L. Ristori, 27 W. J. Robertson, 6 T. Rodrigo, 27 S. Rolli, 38 L. Rosenson, 19 R. Roser, 13 

T. Saab, 14 W. K. Sakumoto, 30 D. Saltzberg, 4 A. Sanson!, 9 L. Santi, 36 H. Sato, 37 
P. Schlabach, 7 E. E. Schmidt, 7 M. P. Schmidt, 41 A. Scott, 4 A. Scribano, 27 S. Segler, 7 

S. Seidel, 22 Y. Seiya, 37 F. Semeria, 2 T. Shah, 19 M. D. Shapiro, 18 N. M. Shaw, 29 
P. F. Shepard, 28 T. Shibayama, 37 M. Shimojima, 37 M. Shochet, 5 J. Siegrist, 18 A. Sill, 35 
P. Sinervo, 14 P. Singh, 13 K. Sliwa, 38 C. Smith, 15 F. D. Snider, 15 J. Spalding, 7 T. Speer, 10 
P. Sphicas, 19 F. Spinella, 27 M. Spiropulu, 11 L. Spiegel, 7 L. Stanco, 25 J. Steele, 40 
A. Stefanini, 27 R. Strohmer, 7 J. Strologas, 13 F. Strumia, 10 D. Stuart, 7 K. Sumorok, 19 
J. Suzuki, 37 T. Suzuki, 37 T. Takahashi, 24 T. Takano, 24 R. Takashima, 12 K. Takikawa, 37 
M. Tanaka, 37 B. Tannenbaum, 22 F. Tartarelli, 27 W. Taylor, 14 M. Tecchio, 20 P. K. Teng, 33 
Y. Teramoto, 24 K. Terashi, 37 S. Tether, 19 D. Theriot, 7 T. L. Thomas, 22 
R. Thurman-Keup, 1 M. Timko, 38 P. Tipton, 30 A. Titov, 31 S. Tkaczyk, 7 D. Toback, 5 
K. Tollefson, 19 A. Tollestrup, 7 H. Toyoda, 24 W. Trischuk, 14 J. F. de Troconiz, 11 S. Truitt, 20 
J. Tseng, 19 N. Turini, 27 T. Uchida, 37 F. Ukegawa, 26 J. Vails, 32 S. C. van den Brink, 28 

S. Vejcik, III, 20 G. Velev, 27 R. Vidal, 7 R. Vilar, 7 D. Vucinic, 19 R. G. Wagner, 1 
R. L. Wagner, 7 J. Wahl, 5 N. B. Wallace, 27 A. M. Walsh, 32 C. Wang, 6 C. H. Wang, 33 
M. J. Wang 33 A. Warburton 14 T. Watanabe 37 T. Watts 32 R. Webb 34 C. Wei 6 



26 

21 



l. u. vvdgnci, j. warn, in. jj. vvaiiduc, n.. 

M. J. Wang, 33 A. Warburton, 14 T. Watanab 
H. Wenzel, 16 W. C. Wester, III, 7 A. B. Wick 
H. H. Williams, 26 P. Wilson, 5 B. L. Winer, 23 D. Winn, iU D. Wolmski/ u J. Wolmski,^ 
S. Worm, 22 X. Wu, 10 J. Wyss, 27 A. Yagil, 7 W. Yao, 18 K. Yasuoka, 37 G. P. Yeh, 7 P. Yeh, 33 
J. Yoh, 7 C. Yosef, 21 T. Yoshida, 24 I. Yu, 7 A. Zanetti, 36 F. Zetti, 27 and S. Zucchelli 2 



(CDF Collaboration) 



2 



* Argonne National Laboratory, Argonne, Illinois 60439 

2 

Istituto Nazionale di Fisica Nucleare, University of Bologna, 1-40127 Bologna, Italy 
3 

Brandeis University, Waltham, Massachusetts 02254 

4 

University of California at Los Angeles, Los Angeles, California 90024 
5 

University of Chicago, Chicago, Illinois 60637 
® Duke University, Durham, North Carolina 27708 

7 

Fermi National Accelerator Laboratory, Batavia, Illinois 60510 



University of Florida, Gainesville, Florida 32611 
ii Frascati, Istituto Nazionale di Fisica Nucleare, 1 
University of Geneva, CH-1211 Geneva 4, Switzerland 



Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, 1-00044 Frascati, Italy 
10 



25 



^ Harvard University, Cambridge, Massachusetts 02138 
12 

Hiroshima University, Higashi- Hiroshima 724, Japan 
13 

University of Illinois, Urbana, Illinois 61801 

14 

Institute of Particle Physics, McGill University, Montreal H3A 2T8, and University of Toronto, 

Toronto M5S 1A7, Canada 

15 

The Johns Hopkins University, Baltimore, Maryland 21218 

^ Institut fur Experimentelle Kernphysik, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
17 

National Laboratory for High Energy Physics (KEK), Tsukuba, Ibaraki 305, Japan 

18 

Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720 

19 

Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
20 

University of Michigan, Ann Arbor, Michigan 48109 

21 

Michigan State University, East Lansing, Michigan 48824 

22 

University of New Mexico, Albuquerque, New Mexico 87131 
23 

The Ohio State University, Columbus, Ohio 43210 
24 

Osaka City University, Osaka 588, Japan 

Universita di Padova, Istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy 
26 



University of Pennsylvania, Philadelphia, Pennsylvania 19104 
27 

Istituto Nazionale di Fisica Nucleare, University and Scuola Normale Superiore of Pisa, 1-56100 Pisa, Italy 
28 

University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
29 

Purdue University, West Lafayette, Indiana 47907 

30 

University of Rochester, Rochester, New York 14627 

31 

Rockefeller University, New York, New York 10021 

32 

Rutgers University, Piscataway, New Jersey 08855 

33 

Academia Sinica, Taipei, Taiwan 11530, Republic of China 
34 

Texas A&M University, College Station, Texas 77843 
35 

Texas Tech University, Lubbock, Texas 79409 

36 

Istituto Nazionale di Fisica Nucleare, University of Trieste /U dine, Italy 
37 

University of Tsukuba, Tsukuba, Ibaraki 315, Japan 
38 

Tufts University, Medford, Massachusetts 02155 
39 

Waseda University, Tokyo 169, Japan 

40 

University of Wisconsin, Madison, Wisconsin 53706 
41 

Yale University, New Haven, Connecticut 06520 

(17 March 1998) 



3 



Abstract 



We report the observations of the decays B + — > -ip(2S) K + and 5° — ► 
ip(2S) K* (892)° in pp collisions at a center-of-mass energy of 1.8 TeV us- 
ing a 110 pb _1 data sample recorded by the Collider Detector at Fermilab. 
We also reconstruct the decays B + -► J/tf;K + and B° -► J/ip K*(892)° and 
measure the six ratios of branching fractions of these four decays. The rel- 
ative branching-fraction results are shown to be consistent with phenomeno- 
logical factorization calculations of hadronic S-meson decays. We use the 
world-average branching fraction B(B + — > J/ipK + ) to derive B{B + — ► 
Y>(2S)if+) = (0.56 ±0.08 ±0.10) x 10" 3 , £(5° -► Y>(2S) K*(892)°) = (0.92 ± 
0.20±0.16) x 10- 3 , and B(B° -► J/^ if *(892)°) = (1.78±0.14±0.29) x 10~ 3 , 
where the first and second uncertainties are statistical and systematic, respec- 
tively. 

PACS Numbers: 13.25.Hw, 13.87.Fh, 14.40.Nd 
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I. INTRODUCTION 



Studies of the decays of bound states of bottom quarks and light antiquarks have proven 
to be one of the most effective ways to explore the decay dynamics of heavy quark systems. 
The branching fractions (B) of the decays of the two lowest-lying bound states, the B + 
and B° mesons, depend on a blend of effects due to the weak and strong interactions. The 
measurements of semileptonic decays of B mesons, where a charged lepton and its corre- 
sponding neutrino are produced, have proven to be useful in the development of theoretical 
models that relate the semileptonic branching fractions to each other [0. In a similar way, 
the measurements of the fully hadronic decays (modes where the 5-meson decay daughters 
are hadrons) have also been shown to provide tests of the theory of heavy quark decay |||3| . 

In this paper, we report the observations of the decays B + — > ip(2S) K + and B° — > 
ip(2S) i^*(892)° in pp collisions at a center-of-mass energy of 1.8 TeV. We also observe the 
decays B + — > J/tpK + and B° — > J/ip if*(892) and use our data to measure the ratios of 
branching fractions of the B + and B° mesons to these four final states. We compare our 
measurements with theoretical branching-fraction ratio predictions. Throughout this paper, 
references to a specific decay mode imply the charge conjugate mode as well. 

As illustrated in the diagram in Fig. [I], all four of the J/ip and i[)(2S) decay modes 
are color-suppressed Cabibbo-favored decays; they can only occur when the W boson's 
hadronic decay products, themselves a color singlet, combine with the charm antiquark 
from the flavor- changing decay and the light spectator quark to form color-singlet charmo- 
nium and strange mesons, respectively. Strong interaction effects, however, are expected to 
modify the dynamics of these decays. The most successful theoretical treatments of such 
decays employ the factorization hypothesis, where the decay of the B meson is described 
by processes that take place on different time scales: short-distance hard-gluon exchange 
and the weak nonleptonic decay of the b quark, followed by longer-distance strong inter- 
actions between the final-state partons that produce the two daughter mesons. The decay 
amplitude is factorized into a product of hadronic currents that reduces to the charmonium 
decay constant and the matrix element for the B — > K hadronic current, which consists 
of several form factors |IJJ2[]. Measurements of the rates and polarization of these decays 
confront the assumptions that underlie the factorization hypothesis in 5-meson decays and 
the calculations involving hadronic form factors. 

Exclusive hadronic decays of B mesons are difficult to observe due to their relatively small 
branching fractions (typically 10~ 4 to 10 -2 ) and the small exclusive branching fractions of 
the subsequent charm daughter decays. However, the large production cross section for 
bottom quarks in pp collisions (in the range of 2-3 fib for quarks with transverse momentum 
Pt > 6 GeV/c and rapidity \yi\ < 1) has made it possible to identify relatively large 
samples of specific decay modes ||]. The decay modes that involve a charmonium daughter 
meson have proven to be most amenable to study, as the decay of the charmonium state 
involving two energetic muons yields a distinctive signature that can be used to identify the 
candidate events. The Collider Detector at Fermilab (CDF) Collaboration has published a 
number of measurements of the properties of B mesons using final states involving a J/if> 
meson, including production cross sections, masses, lifetimes, polarizations, and branching 
fractions 

The CDF Collaboration measured the branching fractions of B + , B°, and B® mesons 
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using five different decay modes, all identified by requiring a J/tj) — > decay [§]. We 

have now completed a more extensive study, incorporating a factor of four more data and 
focusing on the final states identified by the presence of a decay of the form 

B^ip(2S)K (1) 
I — n~ 
I — > J/ipn + 7r~ 
I — 

where B is a B + (B°) meson and K is a (K*(892)°) meson. The if* (892)° meson is 
observed through its decay to the K + it~ final state. 

We first describe in Sees. [Ij and |H| the experiment and data-collection procedures used 
for this measurement. In Sec. [TV], we discuss the event selection procedure and present 
the observed rates of the various 5-meson decays. The necessary efficiency corrections to 
convert the observed decay rates to branching fractions are discussed in Sec. |V| where we 
also detail the systematic uncertainties associated with the measurement. We present our 
results in Sec. |VIJ and offer our conclusions in Sec. |V1J| . 



II. EXPERIMENTAL APPROACH 



In principle, the number of observed events for the decay mode B + — > ip{2S) K + can be 
decomposed into the form 

iV obs (^(25) K + ) = 2 fedt ■ a( P p -&)•/„■ B(B + - ij>(2S) K+) ■ £ ^ K+ , (2) 

and similar forms can be written for the other decays. Here, JCdt is the time-integrated 
luminosity, a(pp — > b) is the bottom antiquark production cross section, and /„ is the 
probability that the fragmentation of a bottom antiquark will result in a B + meson. In a 
similar way, we define fa to be the probability of a bottom antiquark to hadronize and form 
a B° meson. We refer to these probabilities as fragmentation fractions and include in these 
fractions contributions from decays of heavier bottom hadrons into final states containing a 
B + or B° meson. The expression B(B + — > ip{2S) K + ) represents the branching fraction for 
this 5 + -meson decay mode, and e^^ 2S ^ K+ is the acceptance and efficiency of detecting the 
ip(2S) K + final state. The factor of 2 accounts for the possibility of reconstruction of B~ or 
B + mesons in each event. 

The observation of a certain number of 5-meson decays in a specific mode can be con- 
verted into a branching-fraction measurement using an expression similar to Eq. How- 
ever, uncertainties in the bottom-quark production cross section [|J can be avoided by mea- 
suring ratios of branching fractions between I?-meson decay modes, a procedure that also 
results in the beneficial cancellation of several detection and reconstruction efficiencies and 
their associated uncertainties. For example, a branching-fraction ratio involving only J/ip 
charmonium mesons is measured as 

B(B°^J/xl>K*(892)°) = iV obs (J/V;^(892) )-/ M -e J ^^ + 

B(B+ -> J/i)K+) N ohs (J/^jK+) ■ f d ■ e W^*(892)° . B(K*(892)° -> K+ U 
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where N Q b s (J / ip K* (892)°) and N Q \, S {J /ip K + ) denote the observed event yields and 
e J/i>K* (892) an( j e J/ipK+ re p resen t the detector acceptance and reconstruction efficiencies, 
which have several common factors, such as the J/ip branching fraction and reconstruction 
efficiency, that divide out of the ratio. Ratios of branching fractions are also beneficial in 
theoretical studies of these decays since several common factors divide out of the amplitude 
expressions P1JT0| . 



These ratios of branching fractions can be used to estimate absolute branching fractions 
using world-average values for the denominator of the ratios. This is particularly useful for 
those ratios that involve the most precisely known branching fraction, B(B + — ► J/tpK + ), 
in the denominator. We use the world-average value of B(B + — > J/tpK + ) to estimate the 
absolute branching fractions of the other three decay modes. 



III. DATA COLLECTION 
A. The CDF Detector 

The Collider Detector at Fermilab is a multi-purpose detector designed to study 1.8 TeV 



pp collisions produced by the Fermilab Tevatron collider ||11|| . The detector has a coordinate 
system with the z axis along the proton beam direction, the y axis pointing vertically 
upwards, and the x axis pointing horizontally. The polar angle 9 is defined relative to the 
z axis and is the azimuthal angle. Pseudorapidity is defined as rj = — ln[tan(0/2)]. The 
CDF detector surrounds the interaction region with three charged-particle tracking detectors 
immersed in a 1.4 T solenoidal magnetic field. The tracking system is contained within a 
calorimeter system that measures the energy flow of charged and neutral particles out to 
1 77 1 < 4.2. Charged-particle detectors outside the calorimeter are used to identify muon 
candidates. 

The innermost tracking device is a silicon microstrip detector (SVX) located in the region 
between 2.9 and 7.9 cm in radius from the beam axis. The SVX is surrounded by a set of 
time projection chambers (VTX) that measure charged-particle trajectories to a radius of 
22 cm. An 84 layer drift chamber (CTC) measures the particle trajectories in the region 
between 30 and 132 cm in radius from the beam. This tracking system has high efficiency 
for detecting charged particles with momentum transverse to the beam Pt > 0.35 GeV/c 
and 1 77 1 ^ 1.1. Together, the CTC and SVX measure charged-particle transverse momenta 
with a precision of a Pr ~ [(0.0066P T ) 2 + (0.0009P|) 2 ] 1/2 (with P T in units of GeV/c). 

The muon detection system has four of its layers of planar drift chambers separated 
from the interaction point by approximately five interaction lengths of material. To reduce 
the probability of misidentifying penetrating hadrons as muon candidates in the central 
pseudorapidity region \rj\ < 0.7, four more layers of chambers are located outside the magnet 
return yoke (corresponding to a further three interaction lengths of material at 9 = 90°). 
An additional set of chambers is located in the pseudorapidity interval 0.7 < \tj\ < 1.0 to 
extend the polar acceptance. The muon system is capable of detecting muons with Pt ^ 1.4 
GeV/c in a pseudorapidity interval |?7| < 1.0. A three-level trigger system is used to select 
candidate collisions for subsequent study. These and other elements of the CDF detector 
are described in more detail elsewhere |LT] . 
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B. The Data Set 



The data were collected in two running periods, the first extending for nine months 
starting in August 1992, and the second extending for 18 months starting in January 1994. 
Several modifications were made to the CDF detector during the hiatus between these two 
running periods. The most significant of these were the replacement of the silicon microstrip 
detector, the commissioning of a trigger system with greater selectivity, and improvements to 
the data acquisition system. The most notable difference in running conditions resulted from 
a rise in the average instantaneous luminosity of the Tevatron accelerator. The mean instan- 
taneous luminosities during the two periods were 3.5 x 10 30 cm" 2 s _1 and 8.0 x 10 30 cm~ 2 s~ 1 , 
respectively, and the peak instantaneous luminosity exceeded 2.0 x 10 31 cm~ 2 s . The 
time-integrated luminosity of the data sample for the two running periods is ~20 pb _1 and 
~89 pb -1 , respectively. 

Despite the differences in detector configuration during the two running periods, we were 
able to treat the two sets of data as a single sample with a total time-integrated luminosity 
of ~109 pb -1 . This was achieved through the use of nearly identical event reconstruction 
techniques and consistent calibration procedures for data collected during the two running 
periods. 

IV. EVENT SELECTION 

A. The J/ip and ip(2S) Trigger Requirements 

A common feature of the four S-meson decay modes studied here is the presence of a 
fi + fi~ candidate consistent with that arising from the decay of a charmonium (J/ip or ip(2S)) 
state. 

We used a three-level trigger system to identify collisions producing two or more muon 
candidates. The first-level trigger required that two muon candidates be observed in the 
muon system. The first-level trigger track efficiency in the muon system rose from ~40% at 
Pt = 15 GeV/c to ~93% for muons with P T > 3.0 GeV/c. The second-level trigger required 
the detection of at least one charged track in the CTC using the central fast track processor 
(CFT) ||12|| , which performed a partial reconstruction of all charged tracks with a transverse 
momentum exceeding ~2 GeV/ c. The CFT track was required to match within ~8° in <ft of 
the muon candidate. The CFT efficiency rose from ~40% at a muon Pt ~ 2 GeV/c to ~94% 
for P T > 3 GeV/c. The third-level trigger required that two reconstructed CTC tracks be 
matched with two tracks in the muon chambers and that the invariant mass of the dimuon 
pair be between 2.7 and 4.1 GeV/c 2 . The efficiency of the third-level trigger requirement 
was (97 ± 2)% for J/ifj decays passing the first and second-level triggers. Deviations from 
these nominal trigger efficiencies were observed to occur during data acquisition and were 
taken into account in our study. There are 3.1 x 10 6 dimuon candidate events that passed 
the third-level trigger requirements. 

Since the average energy deposition of a muon passing through the calorimeter system 
into the muon chambers was 1.4 GeV, we required that all muon candidates have Pt > 
1.4 GeV/c; however, more stringent criteria on the Pt of the muon candidates were imposed 
by the three-level trigger system, which effectively placed a Pt threshold of approximately 
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2 GeV/ c on each of the two muon candidates. The momentum selection of muon candidates 
by the trigger enhanced the signal yield without introducing large systematic uncertainties, 
since the trigger efficiency was precisely measured. 

In approximately 75% of our selected events, the two muon candidates that were iden- 
tified as charmonium daughters were also the muon candidates identified by the dimuon 
trigger system. In many of the remaining events, an additional muon candidate in the event 
satisfied the dimuon trigger requirements. We included these "volunteers" in our analysis 
in order to maximize the sensitivity of the data sample. 

B. Primary- Vertex and Charged-Particle Reconstruction 

We first identified the location of the pp interaction vertex or vertices using the observed 
tracks reconstructed in the VTX detector. These tracks, when projected back to the beam 
axis, determined the longitudinal locations of these interactions. Due to the high instanta- 
neous luminosities, the number of reconstructed interaction vertices in a given event follows 
a Poisson distribution with a mean of ~2.5. We chose as the longitudinal position of the 
primary pp collision vertex for the dimuon candidates the interaction vertex that was clos- 
est to either one of the muon candidates' intercepts along the beamline. This provided a 
measurement of the primary vertex position with an accuracy of 0.3 cm along the beam 
direction. 

The transverse position of the primary vertex was most accurately determined by using 
the average beam position through the detector and the longitudinal primary vertex position. 
The average trajectory was stable over the period that a given pp beam was stored in the 
Tevatron. The uncertainty in the transverse position of the primary vertex was dominated 
by the transverse size of the beam, which was 25 /xm in both the x and y directions. 

Candidate K + , and n + trajectories were reconstructed in the CTC and extrapolated 
into the SVX detector to identify hits associated with the given track. We required each 
CTC track candidate to be of high quality by stipulating that a candidate track have a 
minimum number of hits in the CTC. We also required that additional SVX information 
consist of at least two out of a possible four hits (for the earlier of the two running periods, 
at least three hits were required). The CTC tracks were also required to pass through most 
of the active volume of the CTC; this was imposed by demanding that the radius of exit 
from the CTC volume of the charged-track trajectory be at least 110 cm. We also required 
that K + and 7r + candidates have a measured transverse momentum > 0.4 GeV/c in order 
that they be reconstructed reliably. 

C. Event Topology Reconstruction 

We performed an event reconstruction using the muon and charged-track candidates. 
The first step of this procedure was to reconstruct a candidate charmonium decay to the 
dimuon final state. The second step was to associate additional tracks with the charmo- 
nium candidate to form a candidate 5-meson decay. During this process, we made various 
selection requirements, described below, to reduce the combinatorial backgrounds and to 
improve the signal-to-noise ratio for the various decay modes. 
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1. Reconstruction of J ftp and if) (2 S) Decays 



In order to reduce the rate of muon candidates from background sources such as K- 
meson decay- in-flight, we required that each muon candidate observed in the muon chambers 
correspond to a CTC track candidate to within three standard deviations of the multiple- 
scattering and measurement uncertainties in both the transverse and longitudinal planes. 

Backgrounds in the dimuon sample were further reduced by performing a least squares 
fit of the CTC tracks associated with the muon candidates under the constraint that they 
originate from a common vertex. We required that the confidence level of this fit be greater 
than 0.01. The dimuon mass distributions for the J/ip and ip(2S) candidates are shown 
in Fig. [|. We fit these distributions to parametrized signal and background lineshapes to 
determine that our sample consists of (4.39 ±0.01) x 10 5 J /if) decays and (1.31 ±0.04) x 10 4 
ip(2S) decays, where we quote only the statistical uncertainties. Candidate dimuon events 
were fit with the additional constraint that the two-muon mass equal the world-average 
masses |Tj| of the J /if) and if>(2S) mesons, 3.09688 and 3.68600 GeV/c 2 , respectively. We 
required that the confidence level of this fit exceed 0.01. This requirement defined our 
inclusive J /if) and if>(2S) data sets. 

Candidates for the decay ip{2S) — > J/if)7r + 7T~ were identified by combining every J /if) 
candidate identified above with pairs of oppositely-charged track candidates that individu- 
ally had to have Pt > 0.4 GeV/c. The two-particle mass of the two pion candidates was 
required to satisfy 0.35 < M(tt + tt~) < 0.61 GeV/c 2 . The lower limit of this mass range 
was motivated by the known dipion mass distribution for ip(2S) decays [0,^,^,0. The 
upper limit corresponded to the maximum dipion mass allowed for ip(2S) decays. The ef- 
ficiency of this criterion to select ip(2S) decays is demonstrated in Fig. ^, where we plot 
the two-pion mass for observed ip{2S) — > J/if> 7r + n~ decays. The background to the ip(2S) 
decays in this distribution has been removed by performing a sideband subtraction using 
background events in the J/ip7i + 7i^ mass distribution. Also shown are curves representing 
a phenomenological prediction JHJ and a pure phase-space calculation. 

To reduce combinatorial backgrounds in the J/if> tt + tt~ candidate search further, an addi- 
tional least-squares fit was performed constraining all four tracks to emanate from a common 
point, the dimuon mass to be equal to the world-average J /if) mass, and the J/ipn + n~ mass 
to be equal to the world-average if>{2S) mass. We required that the confidence level for this 
fit exceed 0.01. The <//■?/> 7r + 7r~ mass distribution for all the candidates prior to impos- 
ing the ip(2S) mass constraint is shown in Fig. |] and illustrates a narrow if>(2S) signal of 
(3.7 ± 0.1) x 10 3 events above the combinatorial background. 

Our subsequent reconstruction identified K + and 7r + meson candidates that were con- 
sistent with arising from one of the four .B-meson decay modes we considered in this study. 



2. Reconstruction of B -Meson Exclusive Decays 

The reconstruction of the .B-meson decay modes required the use of selection criteria 
that reduced the potentially large combinatorial background. The most effective way of 
reducing these backgrounds was to impose minimum Pt requirements on the candidate 
strange-meson daughters, Pt requirements on the .B-meson candidates, and requirements 
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that the decay topology be consistent with that expected from a S-meson decay. Explicit 
particle identification of K + and ir + mesons was not employed in this analysis. 

These additional requirements imposed on the candidates were also selected to be as 
common as possible over the four B-meson decay channels and both decay modes of the 
'0(25') meson to avoid significant systematic uncertainties arising from the estimation of 
5-meson yields across different kinematic regions. The kinematic selection criteria were 
optimized by maximizing the quantity N s /y/N s + Nj,, where N s was the predicted number 
of signal events based on Monte Carlo calculations and iV b was the number of observed 
background events in the 5-meson mass signal region estimated by performing extrapola- 
tions from the sideband regions. This technique avoided the potential bias that could be 
introduced by choosing selection criteria based on the number of observed candidate signal 
events. 

We required that each K + candidate have Pt > 1.5 GeV/c for the reconstruction of 
B + candidates. For the reconstruction of B° candidates, we required that the K*(892)° 
candidate have Pt > 2.0 GeV/c. In the latter case, the K + and tt" daughters from the 
if* (892)° decay were required to have Pt > 0.4 GeV/c. In the K* (892)° reconstruction, all 
possible charged-particle candidate pairs were considered with both mass assignments. A 
track-pair mass assignment was considered a candidate K* (892)° decay if the two-particle 
mass was within 0.0800 GeV/c 2 of the world-average K* (892)° mass of 0.8961 GeV/c 2 . For a 
small fraction of the track-pair candidates, it was possible that both the K + tt~ and ir + K " 
mass assignments fell within this mass window around the K*(892)° pole, resulting in double 
counting of signal events. This double counting was taken into account using a Monte Carlo 
calculation described below. 

To identify S-meson candidates, a least-squares fit was performed on the charged-particle 
tracks associated with the charmonium and strange-meson candidates, constraining them 
to originate from a common decay point with the charmonium-candidate mass constraints 
described above. In addition, this fit constrained the momentum vector of the 5-meson 
candidate to be parallel to its flight path, defined by the measured production and decay 
points, in the transverse plane and required that the B candidate originate from the primary 
interaction vertex. We required the confidence level of this fit to exceed 0.01. We also 
required that the transverse momentum of the B + candidates be greater than 6.0 GeV/c 
and the transverse momentum of the B° candidates be greater than 9.0 GeV/c. These 
different 5-meson Pt selection criteria constitute the largest difference in the kinematic 
requirements between decay modes involving B + and B° final states and are a result of the 
different levels of combinatorial background in the if)(2S) final states. 

Two additional criteria were imposed to reduce combinatorial backgrounds further. In 
the fragmentation of b quarks into B mesons, the meson typically carries most of the energy 
of the quark created in the hard scattering interaction ||18[| . We exploit this fact to suppress 
backgrounds by defining an isolation variable 



where the sum is over charged particles with momentum vectors Pj, contained within a cone 
in rj-<j) space of radius R = J (A0) 2 + (A77) 2 = 1.0 about an axis defined by the direction of 
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the .B-candidate momentum Pg. Track candidates belonging to the B- meson candidate were 
not included in this sum. In order to avoid including charged particles that resulted from 
interactions in the pp collision not associated with the 5-meson candidate, we made the 
sum over only those charged-particle tracks that passed within 5 cm along the z axis of the 
primary interaction location. We expect S-meson decays to have relatively small values of 
Ib, and therefore imposed the requirement Ib < 0.54, which resulted from the optimization 
procedure. We show in Fig. |5| the distribution of Ib for B + — > J /if) K + candidate decays, 
after statistically subtracting the combinatorial background under the B + decay signal using 
background events in the J / if) K + mass-sideband regions. This illustrates that the S-meson 
decays are efficiently identified by this requirement. 

As a final 5-meson selection requirement, we exploited the relatively long lifetimes of 
B mesons and the excellent secondary vertex resolution of the CDF detector to reject those 
events that have short decay lengths 0. We measured the proper decay length for each 
decay candidate, 

_ P T -x T , , 

ct b = — 52 — m B , (5) 

T 

where mg is the mass of the 5-meson candidate, xt is the flight path measured in the 
transverse plane, and Pt is the I?-meson vector transverse momentum. The ctb resolution 
depended primarily on the number of track candidates possessing SVX hit information. We 
required that ctb > 100 fim. 



D. i?-Meson Signals 

The candidate mass distributions for the B + — > J/if>K + and B + — > if)(2S) K + de- 
cays are shown in Fig. |6|. The mass distributions for the B° — > J /if) fT*(892)° and 
B° — > if (2 S) K* (892)° decays are shown in Fig. [7|. We see signal peaks at the B + and 
B° masses in all six decay modes. The distributions describing the if(2S) final states con- 
stitute the first observations of these modes in pp collisions. To obtain the yield of B + can- 
didates, each distribution was fit to a Gaussian signal lineshape with a linear background 
parametrization using a binned maximum likelihood technique. The fits were performed 
over the 5-meson candidate mass region from 5.15 GeV/c 2 to 5.60 GeV/c 2 . The lower edge 
of this range was chosen to avoid possible biases from incompletely reconstructed 5-meson 
decay modes where one or more decay daughters went undetected. 

For the B° candidates, a single Gaussian lineshape did not accurately describe the signal 
due to ambiguities in the K*(892)° daughter (K + and tt~ mesons) mass assignments in 
approximately 25% of the signal events. To correct for this effect, we used a Monte Carlo 
calculation to determine the lineshape for the correct and incorrect mass assignments and 
found that each of these was accurately described with a Gaussian parametrization centered 
on the nominal B° mass. The width of the mass distribution arising from the wrong mass 
assignment was 4-6 times broader than the width of the mass distribution associated with 
the right mass assignment, although the parametrization differed for the three different B° 
decay modes. We list the ratio of the peak amplitudes of the wrong versus right mass 
combinations and the ratio of the widths of these two Gaussians in Table [p. The signal 
yields, masses, and resolutions of the B + and B° candidates are summarized in Table f|. 



12 



V. ACCEPTANCE AND EFFICIENCY MEASUREMENTS 



We used a Monte Carlo calculation of 6-quark production and I?-meson decay, followed by 
a detailed detector simulation to study and measure the kinematic and geometric acceptances 
for each decay mode. We used data to estimate the remaining efficiencies associated with 
the reconstruction algorithms and the event selection criteria. 

An advantage of measuring ratios of branching fractions for similar decay modes is that 
many of the acceptances and efficiencies cancel in the numerator and denominator; however, 
there are several effects that do not cancel completely and have associated with them system- 
atic uncertainties that differ depending on which channels are being compared. Examples of 
these include the effect of the decay-length requirement and the polarization of the S-meson 
decay to vector- vector final states. We have taken these into account and arrived at separate 
estimates of systematic uncertainties for each of the six ratios that were measured. 



A. Acceptance and Trigger Efficiency Measurements 

The Monte Carlo calculation used a model for 6-quark production based on a next-to- 



leading-order QCD calculation [19|. This calculation employed the Martin- Roberts-Stirling 



set DO (MRSDO) parton distribution functions [20|| to model the kinematics of the initial- 



state partons; a value for the 6-quark mass of = 4.75 GeV/c 2 ; and a renormalization scale 



of fi — fiQ = J ml + fcy, where hr is the momentum of the b quark in the plane transverse 
to the directions of the incoming protons. We generated b quarks in the rapidity interval 
\yb\ < 1.1 and with kx > 5.0 GeV/c. These generator limits were chosen so that there were 
no biases in the S-meson kinematic distributions after the application of the kinematic cuts 
used in this analysis. The b quarks were fragmented into B mesons according to the Peterson 
fragmentation function [|18[ with the parameter % defined to be 0.006 |21 . 



The generated B mesons were decayed into the various final states using two-body decay 
kinematics governed by the world-average masses and lifetimes of the daughter particles [JT3j . 



For the ip(2S) — ► J/ipn + n mode, the decay matrix element used was the Pham et al. 
model |TT| , |22|| , a parametrization of which is shown in Fig. [| 

For the J/if> K*(892)° and if)(2S) K*(892)° decays, which involve two vector mesons in the 
final state, we used the world-average longitudinal polarization for the J /ipK* (892)° decay, 
T L /T = 0.78 ± 0.07 0. The polarization has not been measured in B° ip(2S) AT* (892)° 
decays; therefore, we use the value of T L /T for the B° — > J/ip K*(892)° decay but double 
the uncertainty to ±0.14. 

The resulting generated events were passed through a detailed simulation of the CDF 
detector that took into account decays-in-flight, the geometry of all the subdetector elements, 
the interaction of the charged particles with the material in the detector, the resolution of 
the different tracking elements, and the efficiency of the trigger. 

We used a parametrization of the muon-Pr dependence of our first and second-level 
trigger system to determine the trigger efficiency for the various 5-meson decays in which 
the muons that resulted from the charmonium decay were also identified by the trigger's 
dimuon selection criteria. In the remaining events, which corresponded to approximately 
25% of the selected event sample, an additional muon candidate in the event often satisfied 
the dimuon trigger requirements. This led to a ~5% uncertainty in the topology dependence 
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of the relative trigger efficiency. We included this as a separate systematic uncertainty in 
our branching-fraction measurements. 

The resulting geometric acceptances for the different decay modes are presented in Ta- 
ble m, where we list the fraction of decays expected to contribute to the observed signal 
peak. 



B. Reconstruction Efficiencies 

The efficiencies of the subdetectors and reconstruction algorithms used in this analysis to 
identify charged-particle candidates and reject backgrounds were separately estimated using 
information from the collected data together with Monte Carlo calculations. The use of 
data that properly represented the running conditions experienced at the Tevatron Collider 
was particularly important in cases where these efficiencies were expected to depend on the 
instantaneous luminosity. 



1. Track Reconstruction Efficiency 

Because of the different number of charged particles used in the reconstruction of each 
of the exclusive decay modes in this analysis, we had to estimate the efficiency of the 
track reconstruction algorithms when reconstructing three, four, five, or six-prong decays. 
The large number of charged particles associated with an interaction producing a B meson 
(typically ~40 charged particles in the fiducial region of the CTC) and the large number of 
simultaneous interactions created very high hit occupancies in the innermost layers of the 
CTC, reducing their effectiveness. 

We measured the efficiency of the CTC track-finding algorithm by embedding the wire 
hits from one or two Monte Carlo charged-particle tracks into a set of data events identified 
as having a J/ip candidate. This event sample was selected to be representative of the 
inclusive J/ip and ip{2S) data set, taking into account variations in detector configuration 
and instantaneous luminosity. Our embedding procedure used the hit detection efficiencies 
and resolutions observed in the data to simulate the response of the CTC. We then used the 
standard track reconstruction algorithms to seek a reconstructed track helix that formed a 
match with the embedded particle's trajectory. The track matching criterion was imposed 
on a x 2 variable with five degrees of freedom that accounted for parameter correlations 
and measurement uncertainties by employing in its definition the covariance matrix of the 
reconstructed track helix. For an embedded Monte Carlo track in a data event, \ 2 values 
were computed for all of the reconstructed tracks in the event. A track was deemed to be 
a match if it had the lowest x 2 °f an tracks in the event and if this \ 2 value was <500, a 
highly efficient requirement. 

We measured the efficiency for reconstructing a 7r ± meson as a function of the meson's 



Pt, shown in Fig. g| and as a function of several other variables ||23|| . We found the average 
efficiency for 7r ± mesons to be 0.928 ± 0.003 ± 0.026 for meson P? > 0.4 GeV/c, where the 
quoted uncertainties are statistical and systematic, respectively. Because of the use of the 
decay mode i/j(2S) — > J /ip n + n~ , which introduced two additional charged tracks, we also re- 
peated this procedure by embedding pairs of charged tracks in each real event. We found the 
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efficiency for reconstructing the two daughter pion tracks in the decay ip{2S) — > J/ip7i + tt~ 
to be 0.881 ± 0.005 ± 0.043. The systematic uncertainties in these two measurements are 
largely correlated because they were dominated by the uncertainty in the CTC hit efficien- 
cies used in the track embedding procedure. Since the square of the single-track efficiency 
is substantially smaller than the two-track reconstruction efficiency, we concluded that the 
track-finding efficiencies of several charged particles in a single event were correlated, an 
effect that was taken into account in our subsequent efficiency calculations. 



2. Requirements on the Fits to Decay Topologies 

We measured the efficiencies of the constrained-fit confidence-level criteria for the differ- 
ent decay topologies using the observed J/ip and il>(2S) signal yields before and after making 
the confidence-level requirements discussed in Sec. [IV Q . These comparisons resulted in effi- 
ciencies for the constrained fits to the J/xjj — > /i + and the ip(2S) — > /i + fi~ decay hypothe- 
ses that were equal within uncertainties. The efficiency for satisfying the confidence- level 
requirement on the fit employing a vertex constraint was 0.967 ± 0.003 and the efficiency 
for the fit employing the additional charmonium mass constraint was 0.963 ± 0.002. Since 
these efficiencies do not depend on the charmonium parent, they cancel in the measurement 
of the ratios of branching fractions. The uncertainty in these measurements arises from the 
finite statistics of the two charmonium samples. 

The third charmonium topology, i[)(2S) — > J/ip 7r + ix~ , includes two low-momentum 
charged pions in the decay vertex in addition to the muons. We measured the extra in- 
efficiency caused by the vertex-constrained fit to this decay by measuring the loss of ip{2S) 
signal events when making the confidence- level requirement. This resulted in a correction 
factor of 0.834 ± 0.039 for the efficiency of the requirement on the vertex fit alone and a 
factor of 0.945 ± 0.031 corresponding to the additional mass-constrained fit requirement. 



3. Decay-Length Requirement 



The efficiencies of the 5-meson decay-length requirement on the different decay modes 



discussed in Sec. [IV C| were measured by modeling the decay-length resolution, which was 
determined by the number of tracks reconstructed in the SVX detector, and using the world- 
average lifetimes for the B + and B° mesons [13|. We determined the decay-length resolutions 



for events with different numbers of SVX tracks using 5-meson candidate events outside of 
the signal mass region. We then used the Monte Carlo calculation and detector simulation 
to estimate the expected frequencies of events with different numbers of SVX tracks. These 
were then convolved together with the expected exponential decay-length distributions for 
each I?-meson species to obtain a prediction for the observed decay-length distributions. 

The efficiencies of the decay- length requirement were ~0.75 and varied only slightly 
between decay modes. The uncertainties in these efficiencies were dominated by uncertainties 



in the world-average 5-meson lifetimes []T3| and in the number of candidate daughter tracks 
with SVX hit information. 
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4- Isolation Requirement 



The I?-meson isolation criterion discussed in Sec. [IV C| rejected combinatorial back- 
ground. The efficiency of this selection requirement is 0.928 ± 0.054, which we measured by 
estimating the loss of 5-meson decays when applying this criterion for each decay mode. 
We note, for example, that we improve the signal-to-noise ratio in the J/ipK + channel by a 
factor of three when applying this criterion. 

We assume that the isolation criterion is independent of the final-state topology. We 
do not expect this isolation requirement to depend on the specific 5-meson final state, 
since it was defined using only the 5-meson momentum and the momentum of the charged 
particles produced in the production and fragmentation of the b quark. We verified this by 
measuring the efficiencies independently in each channel, but the statistical power of this 
check is limited. 



C. Relative Acceptance and Efficiency Corrections 

For the purpose of determining relative branching fractions, the geometric and kinematic 
acceptance results listed in Table |T| were combined into twelve acceptance ratios involving 
the six reconstructed channels. Table [TV] provides a summary of these relative acceptances 
and their systematic uncertainties. The uncertainties include contributions from the follow- 
ing sources: finite Monte Carlo statistics, variations in the Monte Carlo Pt spectrum with 
changes in the renormalization scale and 6-quark mass, modeling effects in the simulation of 
the second-level trigger, variations in the assumed longitudinal polarization fractions for the 
decays with vector-vector final states, and the uncertainty in the CDF detector simulation. 

Table |V] is similarly structured to present the ratios of the products of the remaining 
efficiencies described in Sec. |V Jjj , namely those associated with track reconstruction, the 
constrained- fit confidence- level criteria, and the proper decay- length requirement. 



VI. RESULTS 

We present our results as six ratios of acceptance-corrected rates of 5-meson decays into 
the four channels. In ratios involving ip(2S) mesons, the results for the two ip(2S) decay 
modes were combined, taking into account the ip{2S) branching fractions. The observed 
numbers of signal events for each decay, listed in Table |J, were divided by the acceptance 
and reconstruction efficiencies listed in Tables [TV] and [V], respectively. The event rates were 



also corrected for the daughter-meson branching fractions listed in Table |V|. Additional sys- 
tematic uncertainties of 4% and 2% were applied to the branching-fraction ratio calculations 
to account for uncertainties in the fitting technique used to estimate event yields and the 
lack of cancellation of the efficiencies of the -B-candidate constrained-fit criteria, respectively. 
When we form the ratios of acceptance-corrected event rates, the 6-quark production cross 
section, time-integrated luminosity, and common efficiencies divide out of the calculations. 
For both the B + and i?°-meson cases, we verified that the event rates for the two different 
ip(2S) decay modes were consistent after taking into account the differences in acceptance 
and reconstruction efficiencies. 
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The measured quantities are the ratios of the product of 6-quark fragmentation fractions 
and the S-meson branching fractions into the specific final state. Thus, our measurements 
can be written as shown in Table |V11| , where the first uncertainty is statistical and the 
second uncertainty is the sum in quadrature of uncertainties in acceptance, efficiency, and 
the charmonium daughter branching fractions. This convention will be employed below 



unless otherwise noted. The ratio in Table [VI I| involving exclusively the J /if) decay mode 



has been measured previously with a subset of these data ||. The present measurement 
supersedes it. 

The results show that the rates of I?-meson decays to the two if)(2S) final states are 
approximately 50% of the rates of the analogous decays to the J/ip final states. In addition, 
we note that the relative rates of vector-vector decays to vector-pseudoscalar decays for the 
if)(2S) and J /if) modes are the same to within relatively large statistical uncertainties. If we 
make the assumption f|[l3| of equal fragmentation fractions, /„ = fa, the vector-vector to 
vector-pseudoscalar decay rate for the J /if) final states, 

B(B° -> J / if) K*( 892)°) 

is now the most precise single measurement of this quantity and provides a constraint on 
theoretical calculations that attempt to describe both this branching-fraction ratio and the 
longitudinal polarization of the vector- vector decay J^|9||LU|] . The analogous result for ip{2S) 
final states may be similarly obtained from Table [VII . 



A. Comparison with Phenomenological Models 

We have compared the results of our measurements to two phenomenological predictions 
for these ratios based on the factorization hypothesis 0,0]. In this class of models, the 
ratios of branching fractions of J /if) and if){2S) decay modes depend largely on the strong 
interaction effects in the final state of each decay, which are modeled by sets of form factors 
specific to each decay. We have chosen two specific models P^f25| , as predictions for the 



branching fractions for ip(2S) and J /if) decays have been made using them. 

The calculation by Neubert et al. [HJ] employs a set of form factors determined using a 
relativistic harmonic oscillator potential model for the meson wave functions and a dipole 
q 2 dependence for most of the form factors, where q 2 is the square of the four-momentum 
exchanged between the B meson and the daughter K meson. (The generic multipole form- 
factor formula is F(q 2 ) = F(0)/(1 — q 2 /m 2 ) n , where m is the pole mass and n = 1 or 2 
for a monopole or dipole dependence, respectively.) The calculation by Deandrea et al. f2"5|| 
determines the form factors empirically by extracting them from semileptonic D-meson 
decays and assuming that they have a monopole q 2 dependence. 

The measured ratios of branching fractions are compared with these two theoretical 
calculations in Fig. |9|, where we have assumed f u = fd- Both models are in agreement with 
the data. However, the results of Neubert et al. are favored overall by the data, whereas the 
calculation of Deandrea et al. predicts a lower branching fraction for the B + — > if)(2S) K + 
decay. 
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B. Absolute Branching Fractions 



We can use the world-average branching fraction [T3 



B(B + -> J /if) K + ) = (1.01 ± 0.14) x 10~ 3 (7) 

and the assumption f u = fa to convert our measurements of branching-fraction ratios into 
absolute branching-fraction measurements for the three other decay modes. There are no 
correlations between our data and the world-average value in Eq. (0), making the determi- 
nation of the resulting uncertainties straightforward. We note that the world-average value 
is based on the assumption of equal f u and fa fragmentation fractions of b quarks produced 
in the decay of the T(4f>) meson. 

The derived branching fractions are therefore 

B(B + -> tp(2S) K + ) = (0.56 ± 0.08 ± 0.10) x KT 3 (8) 
B(B° -> ip(2S) iT(892)°) = (0.92 ± 0.20 ± 0.16) x 10~ 3 (9) 
B(B° -> J/ip #*(892)°) = (1.78 ± 0.14 ± 0.29) x HT 3 . (10) 

The first and second uncertainties are statistical and systematic, respectively. 

The branching fractions for the if)(2S) decay modes are compared with previous mea- 
surements by the ARGUS CLEO and CLEO II J28| Collaborations in Fig. |TG. 



The CDF results are in agreement with these previous measurements and have uncertainties 
that are approximately three times smaller. They represent a significant improvement in 
the knowledge of 5-meson decays to ip{2S) final states. The J /if) A"*(892)° branching frac- 
tion is also in good agreement with the world-average value 0] and a recent measurement 



published by the CLEO II Collaboration pj 



VII. CONCLUSION 

We have made measurements of the rates of the exclusive decays 

B + -> i){2S) K + 

B° -> tl){2S)K*{892)° (11) 

B + -> J/i)K + 

B° -> J/^K*(892)°. 

Using the observed event rates, correcting for the relative efficiencies for these decay 
modes, and assuming equal fragmentation probabilities for B + and -B mesons, we measure 
the following ratios of branching fractions: 



B(B° -> J/i)K*(892f) 
B(B+ J/4>K+) 

B(B+ ^i)(2S) K+) 
B(B+ J/i)K+) 



1.76 ±0.14 ±0.15 



0.558 ±0.082 ±0.056 



18 



B(B°^i/j(2S) if*(892)°) 
B(B+ J/ij)K+) 

B(B+ ^jj(2S) K+) 
B(B° -> J/^A"*(892)°) 

£(E° -> ^(25) i^*(892)°) 
-»• J/^iT* (892)°) 

ff(ff° -> jj{2S) K*(892)°) 
B(B+ ^(2S) #+) 



0.908 ±0.194 ±0.100 (12) 



0.317 ±0.049 ±0.036 



0.515 ±0.113 ±0.052 



1.62 ±0.41 ±0.19. 



These ratios have been compared to two phenomenological calculations that use the factor- 
ization hypothesis. The calculations reproduce the overall features of the trends observed 
in the data. 



We have also used the world-average |13j branching fraction B(B + — > J/ipK + ) to deter- 
mine the absolute branching fractions 

B(B + -> ^(25) = (0.56 ± 0.08 ± 0.10) x 10" 3 
B(B° -> V(25) iT(892)°) = (0.92 ± 0.20 ± 0.16) x 10~ 3 (13) 
B(B° -> J /if>K* (892)°) = (1.78 ±0.14 ±0.29) x 10~ 3 . 

These measurements represent a significant improvement in the knowledge of 5-meson decay 
rates into ip(2S) final states and contribute to an effective test of contemporary I?-meson 
decay models. 



VIII. ACKNOWLEDGMENTS 

We thank the Fermilab staff and the technical staff at the participating institutions for 
their essential contributions to this research. This work was supported by the U. S. Depart- 
ment of Energy and National Science Foundation; the Natural Sciences and Engineering 
Research Council of Canada; the Istituto Nazionale di Fisica Nucleare of Italy; the Ministry 
of Education, Science, and Culture of Japan; the National Science Council of the Republic of 
China; the A. P. Sloan Foundation; the Swiss National Science Foundation; and the German 
Bundesministerium fur Bildung, Wissenschaft, Forschung, und Technologie. 



19 



REFERENCES 



[1] M. Wirbel, B. Stech, and M. Bauer, Z. Phys. C 29, 637 (1985). 
[2] M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C 34, 103 (1987). 

[3] T. E. Browder, K. Honscheid, and D. Pedrini, Annu. Rev. Nucl. Part. Sci. 46, 395 (1996). 
[4] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 75, 1451 (1995). 
[5] F. Abe et al. (CDF Collaboration), Phys. Rev. D 53, 3496 (1996). 
[6] F. Abe et al. (CDF Collaboration), Phys. Rev. D 57, 5382 (1998). 
[7] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 75, 3068 (1995). 
[8] F. Abe et al. (CDF Collaboration), Phys. Rev. D 54, 6596 (1996); Phys. Rev. Lett. 76, 
2015 (1996). 

[9] M. Gourdin, A. N. Kamal, and X. Y. Pham, Phys. Rev. Lett. 73, 3355 (1994); R. Aleksan, 
A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, Phys. Rev. D 51, 6235 (1995). 

[10] A. N. Kamal and A. B. Santra, Phys. Rev. D 51, 1415 (1995); Z. Phys. C 72, 91 (1996). 

[11] F. Abe et al. (CDF Collaboration), Nucl. Instrum. Methods Phys. Res. A 271, 387 (1988); 
D. Amidei et al, Nucl. Instrum. Methods Phys. Res. A 350, 73 (1994); F. Abe et al. 
(CDF Collaboration), Phys. Rev. D 52, 4784 (1995); P. Azzi et al, Nucl. Instrum. Meth- 
ods Phys. Res. A 360, 137 (1995). 

[12] G. W. Foster, J. Freeman, C. Newman-Holmes, and J. Patrick, Nucl. Instrum. Meth- 
ods Phys. Res. A 269, 93 (1988). 

[13] R. M. Barnett et al. (Particle Data Group), Phys. Rev. D 54, 1 (1996) and 1997 off-year 
partial update for the 1998 edition (URL |http: / /pdg. lbl.gov]) . 

[14] G. S. Abrams et al. (Mark I Collaboration), Phys. Rev. Lett. 34, 1181 (1975). 

[15] T. A. Armstrong et al. (Fermilab E760 Collaboration), Phys. Rev. D 55, 1153 (1997). 

[16] J. Schwinger, K. A. Milton, W.-Y. Tsai, and L. L. DeRaad Jr., Proc. Natl. Acad. Sci. USA 
72, 4216 (1975); L. S. Brown and R. N. Cahn, Phys. Rev. Lett. 35, 1 (1975); R. N. Cahn, 
Phys. Rev. D 12, 3559 (1975). 

[17] T. N. Pham, B. Pire, and T. N. Truong, Phys. Lett. 61B, 183 (1976). 

[18] C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, Phys. Rev. D 27, 105 (1983). 

[19] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B327, 49 (1989). 

[20] A. D. Martin, W. J. Stirling, and R. G. Roberts, Phys. Rev. D 47, 867 (1993). 

[21] J. Chrin, Z. Phys. C 36, 163 (1987). 

[22] D. Coffman et al. (Mark III Collaboration), Phys. Rev. Lett. 68, 282 (1992); 69, 3689(E) 
(1992). 

[23] A. Warburton, Ph.D. thesis, University of Toronto, 1998 

(URL |http: / / fnalpubs.fnal.gov/ archive / 1997/thesis/ 1- warburton. ps|) . 

[24] M. Neubert, V. Rieckert, B. Stech, and Q. P. Xu, in Heavy Flavours, edited by A. J. 
Buras and M. Lindner (World Scientific, Singapore, 1992), p. 286. 

[25] A. Deandrea, N. Di Bartolomeo, R. Gatto, and G. Nardulli, Phys. Lett. B 318, 549 (1993); 
R. Casalbuoni et al, Phys. Lett. B 299, 139 (1993). 

[26] H. Albrecht et al. (ARGUS Collaboration), Z. Phys. C 48, 543 (1990). 

[27] D. Bortoletto et al. (CLEO Collaboration), Phys. Rev. D 45, 21 (1992). 

[28] M. S. Alam et al. (CLEO II Collaboration), Phys. Rev. D 50, 43 (1994). 

[29] C. P. Jessop et al. (CLEO II Collaboration), Phys. Rev. Lett. 79, 4533 (1997). 



20 



FIGURES 
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FIG. 1. Diagram of the color-suppressed internal-IF emission mechanism for a B meson (here 
either a bu (B + ) or bd (B°) state) decaying to charmonium (J/ip or tp(2S)) and a strange meson 
(K + or K*(892)°). In this process, the u or d quark is assumed to be a "spectator" of the weak 
interaction. 
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FIG. 2. The dimuon mass distributions for the inclusive J/ip and ip(2S) candidate event samples 
are shown in (a) and (b), respectively. 
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FIG. 3. The observed background-subtracted distribution of the dipion mass (points) in decays 
of the form ip(2S) — > J/ip tt + tt~ . The arrow indicates the minimum mass required in the analysis. 
The solid curve represents a phenomenological prediction due to Pham et al. [jl7| and the broken 
curve describes a pure phase-space distribution. 
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FIG. 4. The J/tpTT + ir mass distribution for the ip(2S) candidates. 
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FIG. 5. The distribution of the isolation variable Ib for candidate B + — > J/tf;K + decays. A 
background subtraction has been performed using the sidebands in the J/ip K + mass distribution. 
The arrow represents the cut below which candidates were accepted. 
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FIG. 6. The J/ipK + mass distribution is shown in (a). The ip(2S) K + mass distributions for 
the two modes ip(2S) — ► /i + /U - and ip(2S) — > J/ipir + ir~ are shown in (b) and (c), respectively. 
The curves are the results of a fit described in the text. 
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FIG. 7. The J/ipK*(892)° mass distribution is shown in (a). The ip(2S) K* (892)° mass dis- 
tributions for the two modes tp(2S) — ► /U + and tp(2S) — ► J/?p7r + iT~ are shown in (b) and (c), 
respectively. The curves are the results of a fit described in the text. 



25 




0.45 0.70 0.95 

P T [GeV/c] 



1.20 



1.2 



0.4 



(b) 



0.0 1 — L 
0.20 



J I I I I I I I L 

0.45 0.70 0.95 

P T [GeV/c] 



1.20 



1.2 



1.0 



0.8 



0.6 



0.4 



0.2 



. ^ 
0.20 



(c) 



J L_l I I U I I I L 

0.45 0.70 0.95 

minimum P T [GeV/c] 



1.20 



FIG. 8. The track reconstruction efficiency for ir + mesons and tt~ mesons as a function of the 
meson Pt is shown in (a) and (b), respectively. The dipion track reconstruction efficiency as a 
function of the lower momentum of the two pions is shown in (c) . 
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FIG. 9. Comparison of the measured branching- fraction ratios (filled circles) with theoretical 
predictions by Neubert et al. [^4| (triangles) and Deandrea et al. |^5| (diamonds). The error bars 
on the measured values represent the statistical and systematic uncertainties added in quadrature. 
In ratios involving B + and B° mesons, f u = fd has been assumed. 
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FIG. 10. A comparison of the derived CDF B(B+ ip(2S) K + ) and B(B° -» ^(25) iT (892)°) 
absolute branching fractions with measurements and limits from the ARGUS pq|, CLEO |2^|, and 



CLEO II [ 28 1 experiments. The hatched bars denote 90% confidence-level (CL.) upper limits and 
the error bars represent the statistical, systematic, and branching-fraction uncertainties added in 
quadrature. 
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TABLES 

TABLE I. The ratios of the peak amplitudes and widths of the Gaussian parametrizations 
describing the wrong to right K-tt mass assignments in the if* (892)° reconstruction. 



£>-Meson Decay 


cc Mode 


Amplitude 


Width 






Ratio 


Ratio 


B° -► J/i> if* (892)° 


J/ip -> n+ (j,- 


0.068 


3.6 


B° -> if)(2S) if* (892)° 


ip(2S) -► |U + 


0.046 


5.8 


B° -> i/j(2S) if* (892)° 


Y>(2S) -»• J/^TT+TT- 


0.043 


6.3 



TABLE II. The numbers of observed signal events, the fitted masses, and the signal widths for 
the six different .B-decay reconstructions. The listed uncertainties are statistical only. 



5-Meson Decay cc Mode Event Yield Mass [MeV/c 2 ] Width [MeV/c 2 ] 



B+ 


-► J ftp if + 


J/V- 






856.7 ±38.3 


5278 ± 1 


14.6 ±0.6 


B+ 


-► V(2S) if + 


^(25) 






71.9 ±13.4 


5281 ± 3 


13.1 ±2.2 


B+ 


-»• V(2S) if + 


Y>(2S) 






37.4 ± 7.4 


5279 ± 2 


11.0 ±2.0 


B° 


-► J/V> if* (892)° 


J/iP- 






378.8 ± 24.8 


5280 ± 1 


13.9 ± 1.1 


B° 


-»• V(25) if* (892)° 


^(25) 




11+11- 


20.9 ± 7.3 


5275 ± 3 


8.2 ±2.4 


B° 


-»• V(25) if*(892)° 


^(25) 




J/lJj 7T+ 7T~ 


29.1 ± 7.5 


5285 ± 4 


11.9 ±3.3 



TABLE III. A summary of the absolute products of the geometric and kinematic acceptances, 
calculated for each decay mode for B mesons produced with kx > 5.0 GeV/c and |y&| < 1.1. The 
uncertainties given are due to Monte Carlo statistics alone. 



S-Meson Decay cc Mode Acceptance (x 10 3 ) 

B + -> Jfip if + J/Y> n+ n~ 19.2 ±0.3 

B+^^{2S)K+ iP(2S)^fi+n- 21.8 ±0.3 

B+^^{2S)K+ iP(2S) -> J/^vr+vr- 6.3 ± 0.2 

B° -► J/iP if* (892)° J/il>-> H + V~ 7.8 ± 0.2 

B° -► V(25) if* (892)° %I){2S) -»• 8.3 ± 0.2 

£° ->■ V(25)if*(892)° ^(25) -»• J/^vt+tt- 3.9 ±0.1 
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TABLE IV. A summary of the geometric and kinematic acceptance ratios with associated 
systematic uncertainties, combined in quadrature, from the following sources: Monte Carlo statis- 
tics, the generated S-meson Pt spectrum, trigger effects, helicity effects, and the CDF detector 
simulation uncertainty. 



Denominator B + ^J/ipK + 
NuuuM'alor > ^ ^ 



B° -► J/%b K* (892)° B + -► ^(25) K+ 

u 



■M M 



B° -► J/^ IT (892)° 



0.405 ± 0.025 



5+ -► V(25) if+ 



1.137 ±0.084 



2.810 ±0.285 



B + -► V(25) 



J/^ 7T + 7T 



0.329 ±0.017 



0.813 ±0.058 



0.289 ± 0.026 



#° V(25) iT(892)° 



0.431 ± 0.044 



1.066 ± 0.092 



0.379 ± 0.030 



B° -► V(25) iT (892)° 

u 



0.201 ±0.017 



0.496 ± 0.032 



0.177 ±0.023 
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TABLE V. A summary of the efficiency-product ratios and their associated systematic un- 
certainties from the following sources: track reconstruction efficiencies, confidence-level criteria 
efficiencies, and decay-length requirement efficiencies. 





Denominator 


5+ J/-0i^ + 


B° -» J/tpK* (892)° 


B + -► i/>(2S) K + 


Numerator 


1 + - 


1 + - 
1 — * i 1 t* 


1 + - 


B°- 


-> J/ip K* (892)° 
1 — >n + (i~ 


0.930 ± 0.030 






B + - 


ip(2S) K + 

1 , 

1 — >P M 


0.980 ± 0.023 


1.05 ±0.03 




B + - 




703 ± 055 


755 ± 062 


717 ± 056 


B°- 


-» V(25) IT (892)° 
' — ► M + M~ 


0.923 ± 0.030 


0.992 ± 0.024 


0.942 ± 0.031 


B°- 


-> V(25) iT (892)° 


0.651 ± 0.053 


0.700 ± 0.055 


0.664 ± 0.054 



TABLE VI. Branching fractions of the daughter- meson decay modes reconstructed in the 
present analysis. The world- average branching fractions were used for the charmonium mesons [13]. 
Note that the branching fraction for the ip(2S) — > fi + fjT decay assumes lepton universality. The 
if* (892)° branching fraction is based on isospin symmetry. 



Decay Mode 


Branching Fraction 


J/ip -> n + n~ 


(6.01 ±0.19) x 10~ 2 


ip{2S) -» n + yr 


(8.5 ± 0.7) x 10~ 3 


ip(2S) -» J/i)K + Tt- 


(3.07 ±0.19) x 10" 1 


K*(892)° -> K+tt- 


2/3 
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TABLE VII. The measured ratios of branching fractions times fragmentation fractions for 
the four decay modes. The systematic uncertainties have been calculated taking into account 
cancellations and correlations in uncertainties. 



Quantity Measured Ratio 

/„ B(B"^J/^K*{892)<-') ]7fi + 01 4 + n] , 

JZ B(B+^J/4>K+) 1./0±U.14±U.15 

B(B+->iP(2S)K+) 
B(B+^J/i>K+) 

fa B(B {> ^(2S)K*(892) a ) 
fu ' B(B+^J/ipK+) 

/„ . B(B+^(2S)K+) n 01 y , n n4q , n no fi 

77 B(B»^J/4>K* (892)0) U.^l/ ±U.U4y ± U.Udb 

%^m^ 0.515 ±0.113 ±0.052 

fe • "gS^iffP 1.02 ±0.41 ±0.10 



0.558 ± 0.082 ± 0.056 
0.908 ±0.194 ±0.100 
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